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ABSTRACT 

Microlensing has proven an effective probe of the structure of the innermost regions of quasars, and an 
important test of accretion disk models. We present light curves of the lensed quasar HE 0435- 1223 in the R 
band and in the ultraviolet, and consider them together with X-ray light curves in two energy bands that are 
presented in a companion paper. Using a Bayesian Monte Carlo method, we constrain the size of the accretion 
disk in the rest-frame near- and far-UV, and constrain for the first time the size of the X-ray emission regions 
in two X-ray energy bands. The /?-band scale size of the accretion disk is about 10 15 53 cm (~46r g ), slightly 
smaller than previous estimates, but larger than is predicted by the standard thin disk model. In the UV, the 
source size is weakly constrained, with a strong prior dependence. The UV to /?-band size ratio is consistent 
with the thin disk model prediction, with large error bars. In soft and hard X-rays, the source size is smaller 
than ^10 15 cm (~13r g ) at 90% confidence. We do not find evidence of structure in the X-ray emission region, 
as the most likely value for the ratio of the hard X-ray size to the soft X-ray size is unity. The simulations do 
not support the idea that quasar flux variability is due to coherent changes in accretion disk area. Finally, we 
find that the most likely value for the mean mass of stars in the lens galaxy is ^0.3M Q , consistent with other 
studies. 

Subject headings: accretion, accretion disks - gravitational lensing: micro - quasars: individual 
(HE 0435 -1223) 



1. INTRODUCTION 

The structure of the innermost regions of quasars and ac- 
tive galactic nuclei, where the X-ray and continuum optical 
and ultraviolet emission originates, is not well understood. 



The thin accretion disk model dShakura & Sunyaev||1973 



|Novikov & Thorne 1973| l is widely used to explain the optica 
and ultraviolet (UV) continuum, but has had difficulty repro - 
ducing its observed spectral shape (e.g., |Blaes et "aL] [2001 ). 
Some progress has been made using quasar variability (e.g., 
|Kelly et al.|2009) |MacLeod et al.|2010l|Dexter & Agol|201 1) 
and "disk reverberation" (e.g., |Collier|2001[ i, but the issue re- 
mains unclear. The structure of the X-ray emitting regions is 
even more mysterious, since the maximum temperature at the 
inner disk edge is too low to produce X-rays. Instead, a jet or 
a hot corona above the disk have been suggested as the X-ray 
source, giving rise to direct and di sk-reflected spectral com - 
ponents (see, e.g., the review by Reynolds & Nowak 2003| >. 
Decisive tests of these models have proven elusive. 

Gravitational microlensing of strongly lensed quasars has 
emerged as an effective technique for measuring the spa- 
tial structure of quasar accretion disks. The stars in the 
foreground lens galaxy magnify (or demagnify) the multiple 
quasar images, causing flux variability that is not correlated 
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as would be expected from intrinsic quasar variability. The 
ratio of the angular size of the emission region to the Einstein 
radius of a star in the lens galaxy determines the strength of 
the variability, with large source sizes smoothing out the light 
curves. Since the accretion disk temperature falls with radius, 
the source is larger at longer wavelengths, and the microlens- 
ing effect is weaker. For a standard thin disk, where the ef- 
fective temperature scales as the (3 = 3/4 power of the radius, 
the characteristic size of the disk varies with wavelength as 
. Several studies have used microlensing to put limits on 
accretion disk sizes, generally finding that they are larger at 
fixed wavelength than predicted by the thin disk model, and 



much larger than sizes calculated us i ng the optical flux ( Poo- 
|ley et al.|2007| |Anguita et al.|2008| |Morgan et al.|2010) . In 
particular, Morgan et al. (2010) find that disk sizes scale with 

2/3 

black hole mass in a manner consistent with the M BH scaling 
expected for a roughly constant Eddington fraction. There are 
also studies of the dependence of the accretion disk size on 
wavelength, which generally find that the size increases with 
wavelength (Poindexter et al.||2008[|Bate et aT|2008[ Eigen 



brod et al. 2008, Floyd et al. 2009; Blackburne et al. 201 
Mosquera et al.||2011). One means of solving the size dis 



crepancy is to have a flatter temperature profile, more like th e 
T cx R 1 ' 2 profile of an irradiated disk (see pvlorgan et al. 2010 1, 
but so far the uncertainties in the slope estimates from these 
studies are too large to distinguish between these models. It 
does appear that scattering on large scales and contamination 
from broad li ne emission have too small an effect to resolve 
the problem dDai et al. 1120 101 Morgan et al. 1120 101 IMosquera 

let aL|201ll >. 

In general, the observed effects of microlensing are stronger 
in X-rays than at optical wavelengths, indicating that the X- 
rays are emitted from a very compact region (e.g., [Pooley 



et al |2007| |Morgan et al.|2008|jChartas et al.|2009[|Dai et al 
2010 1. Recently |Chen et al.l ( |201 lj ) demonstrated that the hard 
X-rays from the quasar Q 2237+0305 are even more strongly 
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microlensed than the soft X-rays, indicating that there may be 
temperature structure within the electron gas that presumably 

generates the X-ray emission. 

At a redshift z s = 1.689 , HE0435-1223 ( |Wisotzki et aT] 
[20021 |Morgan et al. 12005] hereafter HE 0435) is lensed by a 
foreground early-type galaxy at a redshift of zl = 0.46 into 
four images arranged in a nearly perfect cross configuration. 
In this paper, we present light curves of the four images of the 
quasar in the R band and in the UV, and additionally use light 
curves in two X-ray energy bands from a companion paper, 
Chen et al. (2012, in preparation). We use these light curves 
to constrain the size of the source at all four wavelengths. The 
whole is greater than the sum of the parts: the simultaneous 
analysis of all four wavelengths yields better constraints than 
any separate "monochromatic" analysis. In addition, combin- 
ing the UV emission, which is thought to originate in the inner 
regions of the accretion disk, with the optical emission from 
larger radii gives us leverage to probe the change in size of 
the disk with wavelength. Likewise, because we have both 
hard and soft X-ray light curves, we can for the first time put 
a quantitative constraint on the ratio of sizes in the two energy 
bands. Finally, we put weak constraints on the mean mass of 
the stars in the lensing galaxy. We outline the data in Sec- 
tion^ describe our analysis method in Section[3] and discuss 
the results in Section |4] 

2. DATA 

We use eight seasons of optical /?-band photometry result- 
ing from an ongoing monitoring campaign using A Novel 
Double Imaging Camera (ANDICAM, DeP oy et al.|2Q03| > on 
the Small and Moderate Aperture Research Telescope System 
(SMARTS) 1.3-m telescope. This well-sampled light curve is 
supplemented by 10 epochs of UV data from the Wide Field 
Camera 3 (WFC3) aboard the Hubble Space Telescope {HST) 
and 4 epochs of X-ray data from the Chandra X-ray Obser- 
vatory. In the rest frame of the quasar, these wavelengths 
correspond to the near and far UV, and to harder X-rays, 
but for convience we usually refer to them according to their 
observed-frame bands. 

2.1. Optical data 

The magnitudes of the four images of HE 0435 have been 
monitored since 2003 August in the R band (rest frame 
2417 A), with a cadence varying from one day to about a 
week. An automated pipeline subtracts instrumental bias and 
flat-fields the images in the standard way. Three to six im- 
ages are obtained per epoch, and we combine them to in- 
crease the signal to noise ratio and reject cosmic rays. We 
measure the fluxes of the quasar images and the lens galaxy 
using the p oint spread func tion (PSF) fitting method described 
in detail by |Kochanek et al.| ( |2506| . The width of t he PSF is 
measured usin g the SExtractor software package ( |Bertin &| 
|Arnouts||l996| l, and we reject epochs where the full width at 
half maximum (FWHM) is greater than 2 "5 or where the sky 
background is particularly high due to cloud cover or exces- 
sive moonlight. For the handful of remaining epochs where 
the magnitudes look grossly anomalous compared to neigh- 
boring epochs, we re-examine the imaging data and invariably 
find that an uncorrected cosmic ray has collided with one of 
the quasar images. We correct these cases by hand. The light 
curves are given in Table [T] The R band covers a region of 
the redshifted quasar spectrum free of strong emission lines, 
but there may be contamination from Fell emission and the 
B aimer continuum. However, the equivalent width is small 



enough that it should not strongly influence our results (see 
|Daiet al.|2010l >. 

We shift the light curves to correct for the light travel time 
delay between the quasar images. We use delay values A?ba = 
9.5 days, Af B c = 7.2 days, and Af BD = -6.6 days. These val- 
ues are the res ult of fits to the full optical light curves using 
the method of Kochanek et al. ( 2006| l. Though we do not ex- 
haustively explore the range of parameters used for th ese fits, 
the re sulting values are consistent with the values of Courbin 
|et al.| ( |2010[ ), and the intrinsic variability of the quasar is not 
strong enough for small (< 1 day) errors in the time delay cor- 
rections to have noticeable effects on our microlensing analy- 
sis. After being shifted, the light curves of images A, C, and 
D are resampled at the epochs of image B's light curve us- 
ing linear interpolation within each individual season. We do 
not interpolate across seasonal gaps. Finally, we bin each sea- 
son of the light curves to reduce the total number of epochs 
by a factor of 2. This helps to reduce the computational cost 
of our microlensing analysis, and does not significantly de- 
grade our sensitivity to microlensing because its time scales 
are still longer than the sampling of the light curves. The 
delay-corrected flux ratios resulting from this processing are 
shown in Figure [1] Taking the ratio cancels out the intrinsic 
variability of the quasar, leaving only microlensing variabil- 
ity. At this wavelength, the microlensing variability is not dra- 
matic, but can be seen, particularly in the A-B curve between 
the fourth and fifth seasons. 

2.2. UV data 

HE 0435 was observed at 10 epochs between 2009 August 
and 201 1 August using WFC3 aboard the HST. The observa- 
tions used the F275W filter in the UVIS channel, providing 
flux ratios at 2750 A (1023 A in the source rest frame). At 
each epoch there are four exposures of 630 s apiece, for a to- 
tal of 2520 s per epoch. For one epoch, two of the four ex- 
posures were lost to a shutter failure, but we are able to use 
the remaining two exposures. We use the Multidrizzle task 
within the PyRAF software packag^jto combine the four flat- 
fielded images and reject cosmic rays. The world coordinate 
systems (WCSs) provided with the images enables the soft- 
ware to accurately align the images. Examining the bad pixel 
masks produced by this process, we find that several pixels at 
the center of each quasar image are incorrectly marked as cos- 
mic rays, so we manually edit the masks to mark these pixels 
as good, and repeat the final drizzle step. The result is a single 
clean and accurate image for each epoch. At this wavelength, 
the images are remarkably empty, with few sources approach- 
ing the brightness of the lensed quasar. In particular, the lens 
galaxy and quasar host galaxy are not detec ted. 

We use DoPHOT (Sch echter et al.|1993| l to perform aper- 
ture photometry on the quasar components in each epoch, us- 
ing an aperture 0."85 on a side, wide enough to capture es- 
sentially all of the flux in a given quasar component but small 
enough that the other components do not contaminate the pho- 
tometry. We convert the instrumental fluxes into ST magni- 
tudes using the standard header keywords PHOTFLAM and 
PHOTZPT. In this filter, the offset from ST to AB magnitudes 
is otab _ wst = 1-532. The ST magnitudes of the four quasar 
components for all 10 epochs are listed with their formal un- 
certainties in Table [2] and the magnitude differences are plot- 
ted in Figure[T] The magnitudes in Tableware uncorrected for 

6 PyRAF and Multidrizzle are products of the Space Telescope Science 
Institute, which is operated by AURA for NASA. 
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FIG. 1. — Light curves of the flux ratios between the images of HE 0435, expressed in magnitudes. Black circles, magenta diamonds, red crosses, and blue 
crosses show the R-band, UV, soft X-ray, and hard X-ray data, respectively. The curves have been corrected for the lensing time delays as described in the text, 
and the optical light curve has been binned by a factor of two. Any remaining variability is considered to be due to microlensing. The light curves as shown are 
the input to the microlensing simulation. It is clear that the X-ray light curves have significantly more microlensing variability than the optical and UV curves. 



Galactic extinction, which we calculate to be 0.36 magnitudes 
along this sightline, using the Schlegel et al. ( 1998 ) value and 
the |Cardelli et aL] ( |1989| ) R v = 3.1 extinction law. Since we 
are only concerned with magnitude differences, this overall 
offset has no effect on our calculations. At the redshift of this 
quasar, Ly/3 Ovi, and other emission lines overlap with the 
F275W filter, but the passband is wide enough that the con- 
tinuum dominates the flux; we estimate that about 10% of the 
total flux is in the emission lines. Contamination from the 
spatially extended emission line region may bias our results 
toward larg er UV source siz es, but the effect should not be a 
strong one ( |Dai et al.|2010| l. Because the latest epoch of UV 
data extends beyond our optical light curve, and because our 
simulations were begun before it was collected, we use only 
the first 9 epochs of the UV light curves for our microlensing 
simulations. 

We shift the UV light curves of quasar im ages A, C, and D 
according to the time delays given in Section [2~T[ and use lin- 
ear interpolation to resample them contemporaneously with 
B. The time delays between the quasar images are small com- 
pared to the typical separation between epochs, so the adjust- 
ments are small. But there is potential for systematic error be- 
cause of intrinsic quasar variability that is not well-described 
by our linear interpolation, so we broaden the uncertainties 
on the light curves to account for this. To estimate this ex- 
tra uncertainty, we scale and shift the optical light curves so 



that they match as closely as possible the UV light curves, 
and treat them as an estimate of the UV variability. The typ- 
ical absolute difference between the interpolated UV magni- 
tudes and the scaled and shifted optical light curve is ^0.036 
mags, and we add this in quadrature to the formal uncertain- 
ties. Since the statistical uncertainties are so small, this sys- 
tematic contribution dominates the error budget for the UV 
light curves. 

The lack of flux from the lens galaxy in the UV affords us 
an opportunity to investigate the existence of a fifth image of 
the quasar. Lensing theory genetically predicts an odd num- 
ber of quasar images, but in practice lensed quasars almost 
always have either 2 or 4 observed images. The "missing im- 
age" is located at the maximum of the light travel time near 
the center of the lens galaxy, and is strongly dem agnified due 
to the sharp curvature of the potentia l (Rusin & Ma 200f| |Mao] 
|et al. |2001 1 |Keeton|2003") |Rusin et al.|2005| ). This makes cen- 
tral images an interesting probe of both the innermost density 
profiles and central black holes of distant galaxies, but so far 
only two central images have been u nambiguously detected 
dWinn et al.||2004l |Inada et al.||2005| [2008] >. To check for a 
central image at UV wavelengths, we create a deep image of 
HE 0435 by stacking the frames from every epoch using Mul- 
tidrizzle. There are small (~1") offsets between the WCSs 
of each epoch, so we measure the pixel shifts between the 
single-epoch images using cross-correlation and provide them 
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TABLE 1 
R-band Light Curves 



HJD- 2450000 



B 



D 



2863.883 
2871.829 
2877.843 
2884.784 
2891.827 
2899.843 
2906.848 
2916.803 
2919.828 
2926.780 



2.065 ±0.012 
2.140 ± 0.014 
2.134±0.012 
2.141±0.012 
2.172±0.015 
2.171±0.012 
2. 184 ±0.008 
2.192±0.010 
2.177±0.013 
2.186±0.012 



2.632± 0.019 
2.628 ±0.022 
2.644 ±0.017 
2.653 ±0.017 
2.689 ±0.023 
2.728 ±0.018 
2.729±0.011 
2.720± 0.015 
2.759 ±0.020 
2.736± 0.018 



2.592± 0.018 
2.561 ±0.019 
2.670± 0.016 
2.661 ±0.017 
2.685 ±0.022 
2.712±0.017 
2.745 ±0.011 
2.713±0.014 
2.738 ±0.019 
2.643 ±0.016 



2.821 ±0.022 
2.790 ±0.023 
2.803 ±0.018 
2.870± 0.019 
2.831 ±0.024 
2.856± 0.019 
2.891 ±0.012 
2.953 ±0.016 
2.910± 0.021 
2.885 ±0.019 



NOTE. — Light curves are in uncalibrated magnitudes. This table is published in its 
entirety online. A portion is shown here for guidance regarding its form and content. 



TABLE 2 
Ultraviolet Light Curves 



HJD -2450000 






A 


B 


C 


D 


5072.079 


18 


.566 


±0.002 


18.917 ±0.003 


18.867 ±0.003 


19.005 


± 0.003 


5135.646 


18 


.664 


±0.003 


19.004 ±0.003 


18.983 ±0.003 


19.096 


±0.003 


5179.185 


18 


.830 


±0.004 


19.198±0.005 


19. 146 ±0.005 


19.288 


±0.006 


5298.162 


18 


.780 


±0.003 


19. 197 ±0.003 


19.090 ±0.003 


19.247 


±0.004 


5400.013 


18 


.684 


±0.003 


19.064 ±0.003 


18.953 ±0.003 


19.136 


±0.003 


5453.678 


18 


.889 


±0.003 


19.244 ±0.004 


19.204 ±0.003 


19.408 


±0.004 


5550.662 


18 


.992 


±0.003 


19.433 ±0.004 


19.219 ±0.003 


19.676 


±0.004 


5615.296 


18 


.878 


±0.003 


19.323 ±0.004 


19. 104 ±0.003 


19.443 


±0.004 


5662.224 


18 


.742 


±0.003 


19.066 ±0.003 


18.949 ±0.003 


19.119 


±0.003 


5783.750 


18 


.450 


±0.003 


18. 854± 0.003 


18.642 ±0.003 


19.012 


±0.004 



NOTE. — Light curves are in ST magnitudes. 

to Multidrizzle. We use PSF fitting to carefully subtract the 
wings of the four bright quasar images and test for a central 
image. The PSF model is constructed from the (exposure- 
weighted) sum of oversampled and appropriately rotated Tiny 
Tim PSFs for the set of 10 epochs. We fix the relative po- 
sitions of the five quasar images, using the positions of the 
bright four and of the lens galaxy from the CASTLES database 
( |Falco et al.|2001| l, and allow the overall position offset and 
the five normalizations to vary. We rule out a central image 
brighter than 26.3 ST magnitudes at 95% confidence; this im- 
plies a flux ratio relative to the faintest image fs/fo < 0.002. 
This non-detection probably arises from either a massive cen- 
tral black hole or a steep density profile in the lens galaxy, but 
it is difficult to say anything more quantitative from these UV 
observations because of the possibility of extinction within 
the lens galaxy. Likewise, if we had detected flux in the cen- 
tral regions, distinguishing quasar emission from nuclear star 
formation or weak AGN activity in the lens galaxy would have 
been difficult with these observations. 



2.3. X-ray data 

HE 0435 was observed on 2009 /12/07, 2010/07/0 4 
2010/10/29 using C handra/ACIS (| Weisskopf et""al 
20031 



and 



Garmire et al 
each observation 



2002 
ks for 



with an exposure time of 12.}: 
Including an archival o bservation made in 
2006 ( |Blackburne "eTaLl[20lT] |Pooley et ap201T) , we have 
accumulated four Chandra observations of the quasar. The 
details of the reduction of these observation are presented in 
a companion paper (Chen et al. 2012, in preparation). We 
separate the photon events into soft (0.4- 1.3 keV) and hard 
(1.3-8.0keV) energy bands and perform imaging and spec- 
tral analyses of the four images of HE 0435 in these bands 



and in the full energy range. In the spectral analysis, we fit 
the spectra using a power law plus Gaussian emission line 
model, representing the X-ray continuum and the Fe Ka emis- 
sion line, both modified by Galactic absorption and absorption 
at the lens redshift. In the imaging analysis, we use PSF fit- 
ting to determine the photon flux for each image in each epoch 
and each band. After obtaining the best fits from the imag- 
ing analysis, we correct for absorption to obtain the intrinsic 
photon fluxes for each image. The resulting four-epoch X- 
ray light curves for each quasar image in three energy bands 
are presented by Chen et al. (2012), and we use the soft and 
hard band light curves in our simulations. The flux ratios, 
expressed as magnitude differences, are shown in Figure[T] 

Since the X-ray light curves are even more sparsely sam- 
pled than the UV light curves, it seems inappropriate to use 
interpolation to attempt a correction for the systematic errors 
caused by lensing time delays combined with intrinsic quasar 
variability. Instead, we simply leave the X-ray fluxes uncor- 
rected, and broaden the uncertainties to account for this error 
source. Unlike the UV light curves, the statistical uncertain- 
ties are much larger than the systematic errors that we expect 
from time delays. 

3. METHOD 

To constrain the source properties (i.e., the size and shape 
of the quasar emission regions) and lens properties (i.e., the 
mean mass of stars and peculiar velocity), we use many sim- 
ulated light curves. These simulations are compared to the 
observed light curves and evaluated using a chi-square esti- 
mator. Each simulated light curve, which is associated with a 
particular vector in parameter space, is thus assigned a likeli- 
hood. We add them together to create a joint likelihood distri- 
bution, which we then combine with priors on the parameters 
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to generate posterior probability distributions for the quanti- 
ties of interest using Bayes' theorem. This Bayesian Monte 
Carlo approach is laid out in detail by Kochanek ( 200 4]l, and 
we use a version that is updated by Poinde xter & Kocnanek| 
(2010b) to include the random motions of the microlens stars. 

We model the light curves for the four quasar components 
at each of the four wavelengths as an independent intrinsic 
quasar light curve S modified by a time-invariant and achro- 
matic macro-magnification fa due to the overall gravitational 
potential of the lens galaxy, by the magnification due to mi- 
crolensing (5/i,, and by a time-invariant "catch-all" term A/i;. 
This last term represents several possible systematic e ffects, 
such as differential extinction in t he lens galaxy (e.g., |Falco| 
et al.|1 999, Mos quera et al.|201 l|i, millilensing by dark mat 
ubs ' 



ter substructure (e.g., Kochanek & Dalai 2004 Fadely & Kee 



|ton|201 1| , systematic errors in the macro-magnification /z,, or 
low levels of unrecognized contamination of the quasar flux 
by the lens or host galaxy. Expressing these quantities in mag- 
nitudes, this is 

to, (tj,X k ) =S(tj,\ k )+Hi + Afi i (\ k ) + 5[i i (t j7 A(X k )) 

= S(tj,\ k )+nu ot , (1) 

where i enumerates the four quasar components (A, B, C, and 
D), and t j and \ k denote the set of epochs and wavelength in 
our light curves. The function A(\ k ) is the area of the quasar 
emission profile at wavelength A^; it is through this function 
that the microlensing term 5/ij depends on wavelength. 

3.1. Microlensing Simulations 

We simulate the microlensing parts of our model (i.e., the 
S/ii term) using microlensing magnification patt erns created 
using the particle-particle/particle-mes h (P3M, Hockney & 
|Eastwood||l981j ) method detailed by |Kochanek| ( |2004| l. The 
patterns give the microlensing magnification of each quasar 
image as a function of the source's position. The patterns ap- 
pear as a network of intertwined high-magnification caustics 
separated by lagoons of l ower magnification. Like Poindex- 
ter & Kochanek (2010a b), we use dynamic patterns, mean- 
ing that they are recalculated for each epoch. This allows 
us to take into account the random velocities of the stellar 
microlenses as well as the parallactic motion of the earth. 
The patterns are periodic by construction in both the lens and 
source planes, allowing us to wrap light curves and stellar mo- 
tions across the edges, and eliminating edge effects when con- 
volving the patterns. 

The general characteristics of the magnification patterns are 
determined by the shear 7 and the convergence K tot = k s + k*, 
where k s is the convergence due to smoothly distributed mat- 
ter and IS due to the microlens stars. These quantities, as 
well as the macro-magnifications are determined by mod- 
eling the lensin g of the system as a w hole. We use the exten- 
sive models of Koc hanek et al.| ( |2006[ ), who parameterize the 
relative contributions of stars and of the dark matter halo using 
fia i L , the normalization of the stellar component relative to its 
best-fit normalization with no dark matter halo. The particu- 
lar models we use set the effective radius of the stellar com- 
ponent R e to 0"86 and the break radius of the Navarro-Frenk- 
White (NFW) halo r c to 10."0. We use three values of f M / L : 
0.1, 0.3, and 1.0. This allows us to explore a range of values 
for this parameter, and for all of our results we marginalize 
over it. Altho ugh there are re asons to prefer a low value of 
Im/l ( see > e.g., Kochanek et al. 2006b, we conservatively elect 
not to apply a prior on its value. This choice has almost no 



effect on our results. 

We randomly scatter stars across the region near each 
quasar component to create the microlens convergence k*. 
Their masses are drawn from a power-law mass function with 
dN/dM cx M" 1 ' 3 , with a dynamic range of 50 between its 
maximum and minimum mass es. This a reasona ble match to 
stellar mass functions (Poind exter & Ko chanek 2010 b]>, and 
the details should have little effect on our results (e.g.JWyithe 
|et al.|20 00). We use a range of values for the mean mass of the 
stars, with (M/M & ) = 0.03, 0.1, 0.3, 1,3, and 10. The magni- 
fication patterns' outer dimensions are 20 times the Einstein 
radius R E [ n of this mean mass, where 



^Ein - D, 



OS^Ein - 



(AG{M) Dos-Pls 
£>ol 



1/2 



:(5.42x 10 16 cm)(M/M G 



(2) 



and (£>ol,£>os,£> ls ) = (1203, 1746, 1093) Mpc are the angular 
diameter distances from observer to lens, observer to source, 
and lens to source, respectively, for D,m = 0.3, = 0.7, and 
Ho = 70 km s Mpc -1 . The patterns used for the optical wave- 
length simulations are 4096 pixels on a side, so that their 
pixel size is approximately 2.6 x 10 14 (M/M Q ) l l 2 cm. Since 
the mass of the central bl ack hole is ~5x 1O 8 M (based on 
the Civ line width, |Peng et al.|2006[ ), this pixel size is about 
3.6(M/M Q ) 1//2 gravitational radii. Since the emission is more 
compact at UV and X-ray wavelengths, for these simulations 
we create more detailed versions of the same patterns, dou- 
bling the resolution to 8192 pixels. Care must be taken when 
recreating the patterns at higher resolution, because of cer- 
tain details of the pattern creation algorithm. Specifically, the 
shear 7 is adjusted slightly from the input value so that the 
periodic scale is an integer number of pixels in both planes. 
When creating the high-resolution patterns we intentionally 
calculate this adjustment as if the pattern dimension were 
4096 instead of 8192. This results in the desired periodic- 
ity, and avoids small distortions of the new pattern relative to 
the old. 

For each quasar image, the microlensing variability is sim- 
ulated by the movement of the quasar across the pattern, as 
well as the evolution of the pattern itself due to the ran- 
dom motions of the stars. The latter is made possible by 
our dynamic patterns, where each star's position is depen- 
dent on epoch. We use a one-dimensional velocity disper- 
sion of 255 km s" 1 in the lens rest frame, estimated from the 
Einstein radius of a singular isothermal sphere (SIS) model 
of the lens. Treu et al. (2006) have shown this to be a good 
estimator for the stellar velocity dispersion. We assume that 
the lens galaxy has negligible bulk rotational velocity. The 
effective velocity of the source across the pattern is the vec- 
tor sum of the peculiar velocities of the observer, the lens, 
and the source, properly scaled for geometry and cosmolog- 
ical time dilation. W e follow the approach of Poindexter 
& Kochanek| ( |2010b| l, taking for the observer's motion the 
cosmic microwave background (CMB) dipole velocity (from 
Hinshaw et al. 2009) projected perpendicular to the line of 
sight, combined with the parallax due to the earth's orbit, and 
treating the remaining two velocities as normally distributed 
random variables. For HE 0435, the projected CMB veloc- 
ity is (363,-56) km s" 1 East a nd North, respectively. Using 
the empirical interpolation of Mosquera & Kochanek (201 1 ), 
we estimate that the rest-frame peculiar velocity dispersions 
for the lens and the source are ai = 277 and erj = 227 km s" 1 , 
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respectively. Cosmological time dilation and geometric pro- 
jection effects combine to suppress the contribution of the 
source, so that the lens velocity dominates the effective ve- 
locity of the source across the pattern. The dispersion of the 
prior on the effective velocity, projected to the lens plane, is 
given b>[] 



290 km s" 



(3) 



The finite size of the source is simulated by convolving 
the magnification patterns with the source's emission pro- 
file. We use a standard thin disk model that radiates as a 
multi-temperature blackbody, with an effective temperature 
that varies inversely with the j3 = 3/4 power of the radius. We 
neglect the inner edge of the disk. Although it is not guar- 
anteed that this is the correct model for the accretion disk, in 
practice th e details of the source's radial profile are not very 
important (Mortonso n et al.|2005| [Congdon et al. 20071. The 
wavelength-dependent area of the disk A(A) is defined as the 
area within the contour where ksT^s = he /A. We vary this 
area in logarithmic intervals of 0.2 dex, using an adaptive al- 
gorithm to explore the peak of the likelihood function, and 
stopping when the source is too small to be resolved by the 
magnification pattern pixels, when the likelihood falls to a 
very small fraction of the peak val ue, or when 30 values of 
the are a have been explored. Like Poindexter & Kochanek 
(2010a), we also vary the inclination of the disk cos / and the 
position angle of its major axis <j) a . We use ten inclinations, 
with cos / ranging from 1.0 (face-on) to 0.1 (nearly edge-on) 
in steps of 0.1. Similarly, we use 18 major axis position an- 
gles, ranging from 0° to 170°. The inclination and major axis 
position angle do not vary with wavelength, but the area is 
allowed to vary independently in our four bands. 

3.2. Bayesian Monte Carlo 

A handful of parameters determine the exact shape of the 
microlensing light curve: the source parameters (A, cos/, and 
(p a ), the lens parameters ((M) and /m/lX an d a particular start- 
ing point and effective velocity of the path across the magni- 
fica tion pattern. As described in de tail by Kochanek (2004) 
and Poindexter & Kochanek (2010b), we use Monte Carlo in- 
tegration to evaluate the Bayesian integrals over the space of 
these parameters. For each combination of source and lens pa- 
rameters, we evaluate 10 s trial paths across the magnification 
patterns. Each trial is evaluated using the chi-square statistic 



m (tj,Xk) 



2 



<Ji (tj,\k) 



(4) 



where er, are the measurement uncertainties. The intrinsic 
quasar light curve S (tj, A^) and the "catch-all" term Afij(Xt) 
are determined for each trial by minimization of x\- Trials 
with chi-square values above a predefined cutoff value are 
discarded, as their contribution to the likelihood is negligi- 
ble. We do not impose any prior on the intrinsic light curve, 
though in principle one could grant a greater weight to tri- 
als that yield light curves with statistical properties matching 
those of other quasars of similar luminosity and redshift (see, 



7 In their Equation 12, Poindexter & Kochanek ( 2010bl neglect a factor of 
(1 +Zs) in the scaling of as, but this introduces only a negligible error in their 
velocity prior. 



e.g., |Kozlowski et al.|20T0l|MacLeod et al.|201Q] >. We do im- 
pose aGau^sIan - jmor7>^^ mags and a 
variance of (0.1 mag) 2 . It is worth asking whether it is appro- 
priate to allow these quantities to vary independently at vari- 
ous wavelengths, as we have done. It is reasonable to expect 
correlations between wavelengths in the quasar light curve, 
so if we impose a prior on the light curve in future work, it 
ought to include such correlations. Forgoing a prior broadens 
our error bars, and is thus the conservative choice. We make 
a similar conservative choice in allowing A/x, to vary with 
wavelength, since some of its contributors are wavelength- 
dependent (e.g., flux contamination from the lens galaxy at 
optical wavelengths, or differential extinction/absorption by 
the interstellar medium of the lens). 

We are interested in trials that provide simultaneous good 
fits to our optical, UV, and X-ray light curves. In principle, it 
is possible to simultaneously evaluate x 2 f° r a given trial at all 
wavelengths, but given the length of our light curves, particu- 
larly in the optical, we find it less computationally demanding 
to treat the four wavelengths separately. We first amass a col- 
lection of trials that provide reasonable fits to our /?-band data, 
and then resimulate each of these using an independently var- 
ied source area for the UV light curve. This results in a new 
set of trials that fit both the optical and UV data. There is a 
balance between adding and removing trials, since each in- 
put trial generates a number of new trials with a range of UV 
source areas, but some of these are rejected due to the % 2 cut- 
off. Each of the resulting trials has the same parameters as the 
trial from which it was generated, and additionally a value for 
the UV source area and a UV chi-square value and A^of- We 
continue this process for the soft and hard X-ray light curves, 
using the output of one simulation as the input for the next, 
and obtaining in the end a final set of trials with independent 
optical, UV, and soft and hard X-ray source areas. 

We convert chi-square values to likelihoods using anin- 
complete Gamma function, as derived by Kochanek (20041. 
This function is robust to small errors in the light curve un- 
certainties, even with many degrees of freedom; in essence, it 
does not give "extra credit" to solutions with x 2 /Ndoi < 1- We 
adapt it for multiwavelength use, defining the likelihood as 



£ocr 




dof.A- 



1. 



(5) 



0,4 , 



where x\ an d AW,* are the chi-square and degrees of freedom 
at wavelength A^, and f^ k is the factor by which we rescale 
the chi-square to account for errors in the uncertainties. For 
the optical, UV, soft X-ray, and hard X-ray light curves, re- 
spectively, AW,* is 585, 24, 9, and 9, fl k is 2.0, 1.75, 1.25 

and 1.25, and the cutoff values of x\-> a bove which trials are 
discarded, are 4.0, 3.5, 2.5, and 2.5. We multiply the resulting 
likelihood distribution by appropriate priors to generate pos- 
terior probability distributions for the parameters of interest, 
marginalizing over the less-interesting parameters. The prior 
is uniform for the starting position of th e tria l paths, and nor- 
mally distributed as described in Section [3~T| for their effective 
velocities. For the other parameters, we use uniform logarith- 
mic or linear priors. In particular, for the area of the source 
we explore the results of using both priors. 

4. RESULTS AND DISCUSSION 

Figures|2]and[3]show posterior probability distributions for 
the projected area of the quasar. In Figure [2] we plot the joint 
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FIG. 2. — Joint posterior probability distributions for the projected area of the source and the microlens mean mass. The contours show 39%, 86%, and 
99% enclosed probability (note the non-standard values). Black curves denote the logarithmic prior, and gray dashed curves the linear prior. The shaded zone 
indicates the region where the source area is smaller than the pixel size of the magnification patterns. The square of the black hole's gravitational radius (for 
A^BH = 5 X 1O 8 M0) is shown with a vertical line. All other parameters are marginalized. 



distribution of the source area and the mean microlens mass at 
optical and UV wavelengths and in both soft and hard X-ray 
energy bands, which are the result of the combined analysis 
of the light curves at all wavelengths. All other parameters 
have been marginalized. The black contours show the results 
for a logarithmic prior on the source area, and the gray con- 
tours are for a linear prior. Ideally, we would like to have 
results that are independent of the priors, but single-epoch 
and sparsely-sampled microlensing results tend to have some 
prior dependence. In general, a logarithmic prior is appro- 
priate for a positive scale-free parameter like the source size 
(or area), since it makes no sense to give, e.g., the range 1 - 
10 cm ten times the weight of 0.1 - 1 cm. However, on small 
scales between the innermost stable orbit and the size asso- 
ciated with the sampling of the light curve, a linear prior is 
arguably more appropriate. One possible approach would be 
to use a "hybrid" prior such as A /(A + r 2 g ) that acts linearly at 
source areas that are small relative to the square of the black 
hole's gravitational radius and more like a logarithmic prior 
for larger areas. For ease of comparison with other results, 
however, we stay with the traditional uniform priors. With a 
logarithmic prior we cannot put a lower limit on the UV or 
X-ray size, but the distribution increases with decreasing area 
until it hits the resolution limit of our magnification patterns, 



depicted as a shaded region in Figure [2] In the top panels of 
Figure[3]we plot the projection of this joint distribution along 
the vertical axis, marginalizing over (M). The effect of the 
resolution limit can be seen on the logarithmic prior (black) 
curves, causing an artificial decrease in probability at small 
source areas. In the bottom panels of Figure [3] we avoid this 
problem by taking a slice across the joint distribution at fixed 
(M) = O.3M . With a logarithmic prior, we find that in the R 
band log(A/cm 2 ) = 30.87!°;^ (68% confidence), while in the 
soft and hard X-ray bands, respectively, log(A/cm 2 ) < 29.63 
and 29.85 at 90% confidence. Because our resolution limit 
prevents the X-ray distributions from converging, the confi- 
dence level is probably higher than 90%. The UV area has a 
very broad distribution which decreases more or less mono- 
tonically from the resolution limit out to larger areas, with a 
tail extending to larger areas than even the optical distribution. 
Since the linear prior puts a much greater weight on larger 
sources, these distributions are biased high relative to the log- 
arithmic prior, especially in the UV, where the area distribu- 
tion is wide. With this prior, the logarithm of the source areas 
in the R band, the soft X-ray band, and the hard X-ray band are 
31.573$, 29.93^75, and 30.16^, respectively. Since the UV 
distributions are so broad and so prior-dependent, it is proba- 
bly best to view the UV result as inconclusive. In Figures [2] 
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FIG. 3 . — Posterior probability distributions for the projected area of the region producing the observed-frame R band, UV, soft X-ray, and hard X-ray emission. 
The left (right) two panels use a logarithmic (linear) prior. In the top panels, all other variables have been marginalized, including (M) , leading to an artificial 
decrease in probability at small areas (see Figure|2j. In the bottom panels, (M) is fixed at 0.3Mq. In all panels, the square of the black hole's gravitational radius 
is marked with a vertical line. 



and [3] we plot the square of the b lack hole's gravi tational ra- 
dius as a vertical line, using the |Peng et al.| ( |2~006| l mass of 
5 x 10 8 M Q . We can think of this as a hard lower limit on the 
area, keeping in mind that the black hole mass is uncertain by 
a factor of 3 or so. Aside from this we cannot set a lower limit 
on the UV or X-ray size without resorting to the linear prior. 

It is interesting to ask which combinations of wavelengths 
have the greatest power to constrain the area. Obviously, with- 
out the UV light curves we have no information about the cor- 
responding source size, and likewise in X-rays. Nor have we 
simulated either the hard or soft X-ray or UV in the absence 
of the /?-band light curve. But what improvement on the opti- 
cal size is achievable by adding, e.g., just the UV data? Fig- 
ure|4]shows the source area probability distributions obtained 
from each step in our sequence of simulations. The top panels 
show the constraint on the R -band area resulting from model- 
ing only the optical light curve. With a logarithmic prior, it is 
broad enough that it does not converge at small areas. Adding 
additional constraints from the UV (middle panels) and soft 
X-ray (bottom panels) light curves cause the optical size dis- 
tribution to become narrower, particularly in the X-ray case. 
Comparing this figure to the lower panels of Figure [3] it is 
clear that the addition of the X-ray data has significant power 
to exclude small-area solutions for the optical size. This is 



to be expected given the greater variability seen in the X-ray 
band than at other wavelengths (see Figure[TJ. The high X-ray 
variability rules out low-velocity solutions which are the only 
way that a very small optical source can exhibit (relatively) 
small amplitudes of variability. The X-ray and UV data have 
relatively modest effects on each other's distributions. 

The inclination of the accretion disk and the position an- 
gle of its major axis have subtle effects on the microlensing 
light curves, and can only be realistically simulated with dy- 
namic magnification patte rns, as is done for Q 2237 + 0305 by 
Poindexter & Kochanek (2010a). The microlensing variabil- 
ity of HE 0435 is not as prominent as that of Q 2237, and its 
lensing geometry is such that the random stellar motions are 
less important relative to the motion of the source across the 
patterns. For these reasons, it is not surprising that our pos- 
terior probability distributions show no significant preference 
for particular values of the inclination cos i and the major axis 
position angle <f> a . Therefore, we treat these parameters as un- 
determined, and simply marginalize the likelihood over them. 

Since A is its projected area of the accretion disk, it depends 
both on the physical size of the disk and its inclination, in the 
sense that more highly inclined disks must be larger to present 
the same area. When calculating the Bayesian integrals, we 
construct posterior probability functions for the disk's scale 
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FIG. 4. — Posterior probability distributions for the projected area of the source for fixed (M) = 0.3Mq .constrained by various combinations of data. Top 
panels: optical source area, constrained only by the i?-band light curve. Middle panels: optical (solid) and UV (dot-dashed) source areas, constrained by R-band 
and UV light curves. Bottom panels: optical, UV, and soft X-ray (dotted) source areas, constrained by the i?-band, UV, and soft X-ray light curves. The left 
(right) panels use a logarithmic (linear) prior. The square of the black hole's gravitational radius is shown with a vertical line. Compare to the bottom panels of 
FigureB] where the source areas are constrained by all four light curves. 



radius Rs(X), removing the inclination effects. This parame- 
ter is the radius in the thin disk where ksT^r) = hc/X. As- 
suming the accretion disk temperature slope f3 = 3/4 as in the 
thin disk model, it is related to the half-light radius of the ac- 
cretion disk by the factor Ruj/Rs = 2.44. The posterior prob- 
ability distributions for R$ are shown in Figure B] for fixed 
(M) = O.3M . With a logrithmic prior on the scale radii, we 
find an R-band size log(/?s/cm) = 15.53S. By way of com- 
parison, for the adopted black hole mass of 5 x 1O 8 M , an 
Eddington ratio of 0.3, and an accretion efficiency of 0.1, at 
this wavelength thin disk theory predicts a value of 15.12. A 
flux-based esti mate yields a value of 14.74, assuming (3 = 3/4. 
(See, e.g., Mosquera & Kochanek] ( 201 l| l for a description 
of these two estimators.) In the UV we can place a 90% 
confidence upper limit on the size log(/?s/cm) < 15.77, con- 
sistent with both the theory size of 14.63 and the flux size 
of 14.24. The expected theory sizes for the optical and UV 
disk are shown in Figure [5] The soft and hard X-ray radii 
have probability distributions similar to each other, indicat- 
ing sizes smaller than 14.96 and 15.07 at 90% confidence, 
respectively. With the logarithmic prior, the UV and X-ray 
distributions all show a small decline at the smallest radii, but 
this is a resolution artifact similar to the one shown in Fig- 



ure [3] Namely, at the smallest scale radii only face-on disks 
are resolved, so only trials with cos / s» 1 contribute. This 
artifact is only important for the two leftmost data points in 
Figure [5] Our logarithmic -prior estimate for the R band disk 
half-light radius, log<7? 1/2 /cm) = log(2.44fl s /cm) = 15. 92*15, 
is slightly smaller than, but stil l consistent wi t h prev i ous es- 
timates at similar w avel engths. [Morgan et al. (|2010| i, [Black- 
|burne et"aT] ( |201 1[ ) and |Mosquera et al.| ( |201ip respectively 
find values of 16Tt?, 16.04±0.38, and 16.32^. In com- 
paring these results, we note that our present calculations ac- 
count for more sources of systematic uncertainty than these 
other studies, which is why they have comparable uncertin- 
ties despite using superior data. 

Our model parameterizes the dependence of the disk area 
on wavelength using the area at each wavelength, rather than 
a normalizing area and power-law slope. Since the X-rays do 
not originate from thermal thin disk emission, only the optical 
and UV sizes constrain the wavelength slope. Therefore the 
two approaches are equivalent in our case, with two parame- 
ters each. To estimate the dependence of disk area on wave- 
length, we calculate posterior probability distributions for the 
ratio of areas. Figure [6] shows distributions for the logarithm 
of the UV area to optical area ratio, along with similar dis- 
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FIG. 5. — Posterior probability distributions for the scale radius of the 
source at a variety of wavelengths. The top (bottom) panel uses a logarith- 
mic (linear) prior. The arrows mark the observed-frame UV and S-band scale 
radii expected from thin disk theory; the flux estimates are smaller by 0.38 
dex. The vertical line on the left marks the estimated gravitational radius of 
the black hole. The mean mass (M) is fixed at 0.3Mq, and all other parame- 
ters are marginalized. 

tributions for soft X-ray to optical, hard X-ray to optical, and 
hard X-ray to soft X-ray. These distributions use a logarith- 
mic prior only, as a linear prior does not make sense for such 
ratios. We also mark the expected value of the UV to optical 
ratio for the thin disk case where A log A = 8/3 A log A. This 
value is consistent with our probability distribution, which is 
not surprising given the width of the UV area distribution. Al- 
though it is not relevent for the thin disk model, the ratio of 
hard X-ray area to soft X-ray area is interesting because of the 
possibility of structure in the X-ray corona. We do not find ev- 
idence of a difference in sizes: the distribution peaks just at 
unity (though it is broad enough that significant differences in 
the source area at different X-ray energ ies are not rule d out). 
This runs a bit counter to the results of |Chen et al.| ( |201 Tj ), who 
measure X-ray light curves indicating energy-dependent size 
differences in Q 2237 +0305. More densely-sampled X-ray 
light curves would be invaluable in addressing this question. 

Figure [7] shows the posterior probability distribution for the 
mean mass of the stars (M). Like the source area, this mea- 
surement is affected by the resolution limits of our magnifi- 
cation patterns (see Figure [2ji. To avoid this problem, which 
would bias our mass results toward smaller mean masses, we 
use the results of the simulati on that is constrained only by 
the /?-band light curve. Like |Poindexter & Kochan ek (2010b ), 
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FIG. 6. — Posterior probability distributions for the ratios of projected ar- 
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Fig. 7. — Posterior probability distribution for the mean mass (M) of the 
microlens stars. All other parameters are marginalized. This result is taken 
from the optical-only simulation, to avoid biases arising from the lack of 
convergence at small X-ray and UV sizes; see text for details. The solid 
line uses a flat logarithmic prio r, while the dashed line uses the result of 
Poindexter & Kochanek ( 2010b I as a prior. 

we partially break the degeneracy between velocity and mean 
mass by virtue of our dynamic magnification patterns, but as 
we have mentioned, the small random stellar velocities in the 
HE 0435 system reduce the power of this method. Indeed, our 
posterior velocity distribution is indistinguishable from the 
prior, indicating that we have not strongly broken the degener- 
acy. Our mean mass probability distribution is consistent with 
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other e stimates of the mean mass (Poindext er & Kochanek| 
2010b). Since our posterior probability distribution is inde- 



pendent of that of Poindexter & Kochanek (2010b), we mul- 
tiply the two distributions together and plot them as a dashed 
curve in Figure|7] This is equivalent to using their distribution 
as a prior. 



As pointed out by Blackburne & Kochanek (2010), if 
quasar optical variability is due to changes in the effective 
area of the accretion disk, there will be a potentially measur- 
able effect on microlensing light curves that depends on e, the 
power-law slope of the relationship between the source area 
and the intrinsic quasar luminosity A(f) oc L(t) e . In the sim- 
plest model of a multi-temperature blackbody disk where all 
the variability is due to changes in disk size, e = 1 . With our 
dynamic magnification patterns, we have the opportunity to 
test this hypothesis by varying the size of the source accord- 
ing to an estimate of the intrinsic quasar light curve S(t). We 
repeat one of our microlensing simulations, for simplicity fix- 
ing fuiL = 0.3 and (M) = 0.3M Q and restricting ourselves to 
face-on disks. We adjust the source area at each epoch by a 
factor equal to a power e of the source flux, keeping the time- 
averaged source area equal to the fiducial value of the area 
parameter. We allow e to take values ranging from -1.0 to 
2.0 in steps of 0.5. The intrinsic light curve is estimated by 
averaging the source light curves corresponding to the best- 
fitting trials from our main optical simulations. The resulting 
probability distribution for the exponent e allows for the en- 
tire simulated range, but favors the lower values. This result 
argues against the idea that coherent changes in the area of the 
accretion disk account fo r quasar v ariability , and contr adicts 
the preliminary tests of Blackburne & Kochanek (2010 ), who 
find qualitative evidence that the area and luminosity of this 
quasar are positively correlated. 

In future microlensing studies using light curves at sev- 
eral UV, optical, and infrared (IR) wavelengths, it would be 
preferable to parameterize the wavelength dependence of the 



accreti on disk size using a power law, as in |Poindexter et al. 
(2008), rather than using independent size estimates at each 
wavelength, as we do in this work. The difficulty with this ap- 
proach is the necessity of modeling all the light curves simul- 
taneously, a task that can quickly outpace our computational 
(specifically, memory) capacity unless the algorithm is care- 
fully designed. When there are contemporaneous observa- 
tions at multiple wavelengths, we can save some memory by 
reusing magnification patterns. Another promising strategy 
would be to separate the data into a densely-sampled single- 
wavelength light curve and a second multi-wavelength light 
curve constructed from relatively sparse light curves at the re- 
maining wavelengths. This would allow the successful trials 
from the single-wavelength simulation to be used as input for 



a second, multiple-wavelength simulation where the power- 
law slope a! log A/a! log A would be allowed to vary. This is 
functionally equivalent to the approach in the previous para- 
graph, but with the wavelengths in place of the estimated in- 
trinsic light curve. 

Based on eight seasons of 7?-band monitoring data, nine 
epochs of UV photometry from HST, and four epochs of X- 
ray fluxes in two energy bands from Chandra, we have put 
limits on the projected area and scale radius Rs of the accre- 
tion disk of HE 0435 -1223. Using logarithmic priors, we 
find in the observed R band that the scale radius log(R s /cm) 
is 15.53^5, or about 46 gravitational radii assuming a black 
hole mass M B h = 5 x 1O 8 M (based on a Civ line width mea- 



surement; Peng et al. 2006). This large disk size confirms 



previous results (|Morgan et al.|2010| [Blackburne et al.|2011 



Mosquera et al. 2011), being about 0.4 dex larger than the size 



predicted by thin disk theory, and 0.8 dex larger than the size 
estimated from the optical flux. The UV size is more poorly 
constrained, but we can say that it is most likely smaller than 
the optical size. In soft (0.4 - 1 .3 keV) and hard (1.3-8 keV) 
X-ray energy bands respectively, we can set upper limits of 
log(/?s/cm) = 14.96 and 15.07, or 12 and 16 gravitational 
radii, with at least 90% confidence. We also compute prob- 
ability distributions for the area ratios between various wave- 
lengths. With a very broad distribution, the UV/optical ratio 
is consistent with both thin disk theory and with previous esti- 
mates of the wavelength slope, and the hard X-ray/soft X-ray 
area ratio is most consistent with unity. Adopting a range of 
values for the mean mass (M) of the stars in the lens galaxy, 
we find results consistent with earlier studies, with a most 
likely range of 0.1 - 1 M . Finally, a test of the power-law de- 
pendence of the source area on the instantaneous luminosity 
of the quasar does not support the idea that quasar variabil- 
ity is a simple result of coherent changes in the area of the 
accretion disk. 
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